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HYPERSONIC AERODYNAMIC CHARACTERISTICS OF A FAMILY 
OF POWER-LAW, WING-BODY CONFIGURATIONS 

By James C. Townsend 
Langley Research Center 

SUMMARY 

The configurations analyzed are half-axisymmetric, power-law bodies surmounted 
by thin, flat wings. The wing planform matches the body shock wave shape. Analytic 
solutions of the hypersonic small disturbance equations form a basis for calculating the 
longitudinal aerodynamic characteristics. Approximate boundary-layer displacement 
effects on the body and wing upper surface are included. Skin friction is estimated by 
using compressible, laminar boundary-layer solutions. By using an effective body shape, 
the method is extended to small angles of attack. Three basic theoretical assumptions 
are made: (1) the body is slender, (2) the shock wave is strong, and (3) the Mach number 

is large. In comparisons with available experimental data, good agreement was obtained 
when these assumptions were satisfied. The method is also used to estimate the effects 
of power law, fineness ratio, and Mach number variations at full-scale conditions. The 
implementing computer program is included. 

INTRODUCTION 

Much research has been devoted to the hypersonic flow about half bodies of revolution 
mounted beneath a thin wing. Theoretical studies (refs. 1 to 3) and experimental work 
(refs. 4 to 7) show that with half-cone bodies these configurations combine good stability 
characteristics with high values of maximum lift-drag ratio. Replacing the conical 
bodies with those having power-law profiles generates a larger class of configurations 
and one which is more representative of aircraft shapes. Low wave-drag bodies in the 
hypersonic regime are generated by power-law curves with exponent in the range 0.5 to 
0.8. (See refs. 8 to 12.) These bodies have the additional advantage of better volume 
distribution than cones. 

The purpose of this study was to develop a method for calculating the longitudinal 
aerodynamic characteristics of power-law bodies with reflection- plane wings. The method 
applies to configurations consisting of half of an axisymmetric power-law body mounted 
beneath a thin wing whose planform matches the theoretical body shock shape at zero angle 



of attack. Small -disturbance theory, with small perturbations for Mach number and 
boundary-layer displacement effects, provides a means for calculating the pressure field 
and shock-wave shape. This pressure field is integrated analytically to obtain the forces 
and moment on the body. Small angles of attack are simulated, and laminar skin friction 
is calculated. The computer programs which have been written to implement this method 
are presented in an appendix. 

SYMBOLS 


shock-wave perturbation constant 


Chapman- Rubesin constant, 


Mw/M c 

T w / T c 


axial-force coefficient, 


Axial force 


drag coefficient, 

3oo S 


lift coefficient, 

q S 


pitching- moment coefficient, Pairing moment 
normal-force coefficient, Normal force 


normal-force coefficient of body 


normal-force coefficient of wing 


p - p 

pressure coefficient, — = — — 

<lco 


mean aerodynamic chord, taken as 5]-, = 


m + 2 


constant in boundary-layer displacement thickness 
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F 

f 

I 

J 

L 

l 

m 

P 

h 

R 

^ 00 , l 

r 

S 

s c 

T 

Uoo 


similarity static-pressure variable 

fineness parameter, — 

r b,B 

boundary- layer profile parameter 
integral of F from body to shock 
lift 
length 

free- stream Mach number 

exponent of power- law body shape 

dimensionless static pressure, — ^ — 

26 ^ 

average wing upper surface pressure 
free- stream dynamic pressure 

p 

dimensionless shock-wave radius, — 

6 1 

PjUoo 

free- stream Reynolds number, 

Moo 

f 

dimensionless radial coordinate, — 

5Z 

projected planform area 

distance from nose to upper surface center of pressure 
temperature 

free-stream axial velocity 
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‘-ft 


volume of body 


dimensionless axial coordinate, ~ 

* l 


a 


angle of attack relative to body axis 


y 

6 


ratio of specific heats 

Rq D I 

shock-wave slope parameter, 6 = — 2— = 

l % 

dimensionless boundary- layer displacement thickness, — 

6Z 


* 


small perturbation parameter for Mach number, 


(5MJ' 


small perturbation parameter for boundary- layer displacement, 


V 

9 


similarity form of radial coordinate, — 

R 0 

shock-wave angle 


f e 

v - ~ 


viscosity coefficient 


similarity form of axial coordinate 
dimensionless density, — 


Subscripts: 


B at base of configuration, x = 1 

b body 


d6*/d£ 

dr b/d| 
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e effective body shape 

max maximum 

(L/D) max maximum lift- drag ratio 

0 zero- order similarity solution ^ — o) 

1 first-order similarity solution ( e ^ « l) 

oo free- stream value 

w wall 

An asterisk denotes that the quantity includes boundary -layer displacement effect. 

A bar over a symbol denotes a dimensional quantity. 

THEORY 

The method applies to the general configuration shown in figure 1(a). It consists of 
one-half of a body of revolution mounted beneath a thin wing at an angle of attack of 0°. 

By assumption, the wing acts as an endplate to maintain the axial symmetry of the flow 
about the body. The wing planform matches the shock-wave shape about the full body, 
and the body pressure field acts on the wing to provide additional lift. The method is put 
together from a series of pieces in order to arrive at the final aerodynamic coefficients. 
The basis for the development is the result in hypersonic slender-body theory that for 
power-law bodies, there are similarity solutions to the inviscid flow equations in the hyper- 
sonic limit. (See ref. 13.) Independent small perturbations are made to account for Mach 
number effects and for laminar boundary-layer displacements. (See ref. 14.) To simulate 
the effects of small angles of attack, a simple substitution of an effective body is made. 

The resulting equation for the pressure distribution is integrated analytically to obtain the 
pressure forces and moments on the body. Then the laminar skin-friction drag is calcu- 
lated by using the analytic pressure distribution. The development outlined is explained 
in more detail in the following sections. 

Inviscid, Power-Law Body Solution 

If, in inviscid hypersonic flow about a slender body, the velocity changes in the free- 
stream direction are neglected compared with the transverse flow velocities, the hyper- 
sonic small-disturbance equations result. When a strong, power-law shock wave occurs 
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in such a flow at infinite Mach number, these equations indicate that the body generating 
the shock also has a power-law shape. (See ref. 14.) Thus in figure 1(a), the shock wave 

^r- = 6 j^fj is generated by the body ^ . In dimensionless form (fig. 1(b)) these 

relations become Rq = x m and r b = ^ x m . When the dependent variables are expresser 

in terms of the slope of the shock wave, the axial variations may be separated from the 
radial variations of the variables to obtain similarity equations. Thus in similarity varia- 
bles | = x and T) = t/Rq; Rq - | m , r b = 7 / b £ m , and the dimensionless pressure field is 

p 0 = F 0 = m 2 FQ( 7 j)(;^ m-1 ). Here Fq(ii) is found by solving a set of ordinary 
differential equations in 77 . (See refs. 14 and 15.) 

In order to relax the restriction to infinite Mach number, Kubota (ref. 14) applied 
a small perturbation procedure. This procedure results in the following first-order pres- 
sure distribution and shock- wave shape about the power -law body r b = 

Pjd,?]) = m 2 F 0 (7j) $2(m-l) + e 1 m 2 F 1 (rj) (1) 


»,«> = r L 1 * (2) 

Here = (M^) 2 is a small parameter corresponding to the hypersonic strong shock 
assumption M ro sin 0 » 1. The necessity of simultaneously satisfying S 2 « 1 and 
€j 2 « 1 puts a strong requirement on the Mach number. That is, the present method is 
limited to M m » 1 so that with a slender body 5 2 « 1 , the parameter e j is still 
small. Figure 2 shows the relationship between 6 , e lt and M m and can be used to 
check on ej for a given Mach number and body shape (by noting that 6 = l/fTj b ). 

The perturbed pressure variable F^(tj) and the shock- wave displacement constant 
a l are found from a second set of ordinary differential equations. (See refs. 14 and 15.) 
The two sets of differential equations involve only m and y as parameters, and thus 
they can be solved over the needed range of values and the tabulated results used in appli- 
cations to flow problems. Table I and figure 3 present the results needed for the current 
application. They were found by numerical integration techniques similar to those 
described in reference 16. (With acess to modern digital computers, the exact numerical 
computation has become at least as easy to carry out as the approximate technique which 
gives ref. 16 its title.) The integrated pressure J 0 and pressure perturbation Jj are 

defined as Jo = \ F o^) d7 7 and Jl 5 \ F l(v) dv‘> they will be applied to the wing 

. . '% J % 
undersurface. 
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Corrections for Boundary-Layer Displacement 

In order to make a corrected approximation accounting for laminar boundary-layer 
growth, a perturbed body shape r b *(x) = r b + 6* is used. The displacement thickness of 

the boundary layer is given by 5* = 2m?7 k - Ex 3 /2-m (based on a result from ref. 17 for 

3 - 2m 

adiabatic wall conditions). Here E s 2Lli. M^Irv — 3 ~ 2m ■ ■ ■£ S, , > which is 

l/2y 2m 2 V4m - l| R »l r o(%) 

very small for large Reynolds numbers. By using the appropriate value for I (the sum 
of the transformed displacement and momentum thicknesses in refs. 18 and 19), this 
relation for 5* may be applied for any constant wall temperature. If the flow outside 
the boundary layer is considered to be the inviscid flow about the "perturbed body" 
r b *(x), the corresponding pressure distribution and shock shape are approximated as 
follows. In terms of the body radius r^ = 'bb| m > equations (1) and (2) become 


P x = m 2 F 0 ( + €l m 2 F 1 (7/) = + e 1 m 2 F 1 (r?) 


_ t m L 2(l-m)l r b , / 1 dr b 

R, = £ 1 +a 1 e 1 4 = — 1 +aiei - 

1 L 1 1 J % 1 X K d£ 


By replacing r^ by r b * (as in ref. 14), these equations become 


/ *\ 2 

(4,??) = £2^ [ — ) + e 1 m 2 F 1 (7j) = m 2 i 2(m_1) (l + €*) 2 F 0 (7?) + e^Fjfa) 

% z V d W 


r h * / 1 dr. * 

R,*U) » — 1 + a iei(— - 37 - 
1 % 1 l \ m % d £ 




+ terms of order e * 


I -2m 


Here e * = ^ = E£ 2 « 1 except in a small region near the nose when m > j. 

dr b /d£ 4 
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Simulation of Angle of Attack 


The pressure distribution along the pitch plane of the body at angle of attack is 
assumed to be the same as that about an equivalent axisymmetric body. This effective 
body is at zero angle of attack. It has a power -law profile which closely matches the 
windward element in the plane of symmetry of the actual body at angle of attack. Figure 4 
shows the relation between the real and effective bodies, and the following expressions are 
used to obtain the effective body parameters 


x e = x cos a - r^ sin a 


r b e = x sin a + r b cos a 


f = 


1 - -j tan o' 


b,e,B tan a + j 


1 1 + f tan a 




(5) 


(tan a « 1) 


and 


so that 


m e = 


_ 1 °ge( f b,2,e/ f b.l,e) 


log, 




m £ 


•b,e w jJV 
^e \/e j 


Here xj and x 2 31,6 points selected to provide a good approximation. A lower limit 
of m e =0.51 was set to avoid computational problems associated with the theoretical 
limit as m — 0.5. 

The approximating pressure distribution along the body at angle of attack is then 


l,ei^b,e 


m. 


1 + 2i^E e x e 2 


#-2m e 


'°(\e) x e ' ' + «l F l(V) 


( 6 ) 


where 


v = -$■ = 


f 1 + f tan a 


* , (±\- , iV 
\% / 
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E with m e replacing m. |The factor k appears since 
By following reference 6 which uses a similar angle -of-attack 

method for half-cone wing configurations, this pressure distribution is applied over the 
entire body surface. The pressure distribution under the wing from rj^ = 77 = 1 is 
assumed to be the same as that in the flow field of the effective body from 7j b>e = v = 1. 
Since this equivalent body approach does not attempt to account for the actual flow under 
the wing, its use necessarily limits the present method to very small angles of attack. 

For this reason all calculated results presented are in the range ^ = Instead of 

calculating the wing upper surface pressure in detail, an average pressure is used. This 
value is taken from the charts of reference 20, which includes viscous -interaction effects 
on the pressure and skin friction on delta wings at angle of attack in hypersonic flow. 

Since the viscous effects are approximately proportional to x - !/ 2 , delta-wing results for 

which ^yx _1 / 2 drdx and the span equal to those of the power-law wing are used. The 
base pressure is set equal to free-stream static pressure p^. 


Skin Friction 

The skin-friction contribution is the remaining term of the axial-force coefficient to 
be evaluated. In this report laminar boundary layers are assumed for all calculations. 
The wetted area is divided into the body surface, wing upper surface, and the exposed part 
of the wing underside, each of which is treated separately. For the skin-friction calcula- 
tions for the body and the wing lower surface, the longitudinal pressure distribution is 
modified in the nose region by keeping a higher order term in the pressure equation. 

These calculations then use a scheme given in reference 18 for incompressible laminar 
boundary layers. Two transformations of the independent variables allow its use with the 
two-dimensional compressible laminar-boundary-layer similar solutions of reference 19 
for the present cases. For the body, the Mangier transformation (ref. 18) changes the 
axial coordinate to that for an equivalent two-dimensional body. For the relatively small 
exposed-wing undersurface, a simplified flow model is applied, that is, streamlines are 
taken as parallel to the body surface, and the pressure is taken as varying parabolically 
from the body to the shock wave. In both cases the Stewartson transformation (ref. 19) 
changes the surface length and exterior velocity distribution to the form for an equivalent 
incompressible flow; the method of reference 18 is then applied. For the wing upper 
surface, the average skin friction from the appropriate charts of reference 20 is used just 
as for the upper surface pressure. 

Longitudinal Aerodynamic Coefficients 

Integrations of the appropriate components of the surface pressures over the body 
and wing give expressions for the axial force, the normal force on the body and on the 
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wing, and the pitching moment. In coefficient form the expressions are (to first order in 
ej and e*): 


c A = 


-N,b 


jr(m + l)S 2 m p 2 


s b 

2m 2 6 2 ( 


1 4^ 2 Ep . £1K , . 1 

+ m - 1) + 4m - 1 mFo K,e) + ~ ~ 


2(m e + m - 1) 


^TTT + Srf ) (m + 1)F o(%,e) + MXe) - “T 


rm e 


i+ C 


a,f 


_S I \2m 
Sb 

2*2 


Sb 


1 2 


2mg 6‘ 

N " k J- UV 2ra e+*n-l 2m + 1/ V 1 - V b>e 


n a i 


2iru - m + 1 




°> e y + 1 
1 - % 


Y m e 


4mE 


(3 - 2m)(4me - 2m + 1) 


ym. 


C -C + c Pu ' P « 

- L N,b + c N,w n 




(m + l)(m + 2)m e 2 6 2 f j / j 


Sc 

Vb 


4^ 2 E, 


f2 \2nig + 3m - 2 6 m 


~J mF oK, e ) 


3f 2 


F lKe) 


ym 


e J 


1 

1 + 


4^2e, 


2m e + m 2m + 3 


F 0 h bp ) + 


2 

4mE 

(y + t^b 

(3 - 2m)(4m e - 2m + 3) 


e t a l 


ejK 


m + 2 


■K,e) 


1 - % J l,e 1 


1 - 


^b.e ^b ym e \ 


gu ~ ^ rj ° s c 

C 


The pitching- moment reference center is at the nose of the body (x = 0). For the corre- 
sponding zero-order and inviscid relations, set ej = 0 and E = 0, respectively. Note 
that the factors and m + 1 are associated with the actual planform area used in 
normalizing the coefficients. 

These equations have been programed for calculation by a high-speed digital com- 
puter. The program includes the skin-friction calculations on the body and the wing 
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undersurface. There is a subsidiary program which computes the parameters required 
to get the upper surface pressure and skin friction from reference 20. The appendix 
presents an exposition of these programs. The basic program requires only 21 600 octal 
storage locations after compilation on the Control Data Corporation 6600 computer at 
the Langley Research Center and runs an average case with 11 angles of attack in about 
3 seconds of central processor time. 

DISCUSSION OF RESULTS 
Evaluation of Method 

There have been no reported comprehensive experimental evaluations of the power- 
law, wing-body configurations to which the theoretical analysis applies. The data avail- 
able fall into two groups: (1) drag of complete power-law bodies of revolution (no wing) 

at several Mach numbers and fineness ratios, and (2) aerodynamic characteristics of 
conical (m = 1.0) wing-body configurations. Only a small part of these data satisfy the 
high Mach number, the slender body, and the strong shock criteria required for strict 
application of the similarity theory. Data for which the criteria are not well satisfied 
can be used to determine the limits for practical application of the method. 

Power-law bodies of revolution.- The zero-angle-of-attack drag of these bodies is 
already calculated as part of the present method. Figure 5 contains four sets of compari- 
sons with experimental data. The drag coefficients have been based on the length squared 
as reference area in each case to form a uniform basis of comparison. In parts (a) and (b) 
of figure 5, the ratio V/i 3 was held constant and yielded a small variation in the fine- 
ness ratio as the power-law exponent was varied to obtain the different bodies for the tests. 
Figure 5{a) is for tests at Mach 21.6 in helium (ref. 12). The agreement is very good. 

The coefficients in figure 5(b) are for tests at Mach 10.03 in air (ref. 21); the calculations 
are in good agreement with experiment. Figure 5(c) shows good agreement at Mach 10.35 
for a series of power-law bodies having nearly equal fineness ratios. In figure 5(d) the 
data for the same bodies at Mach 5.96 is not predicted. 

The range of agreement obtained in figure 5 should be considered in light of the 
basic assumptions of the theory as discussed in the previous section. For this reason 
the pertinent parameters are shown in the legends of figure 5 and also in figure 2. Since 
6 « 1 for all cases, the slender body condition is well satisfied. The hypersonic assump- 
tion (m » 1 ) is generally considered to be satisfied for M > 5 and so should not cause 

\ CO / / 2 \ 

the discrepancies in figure 5. However, the strong shock assumption « lj is satis- 
fied only for figure 5(a), where the agreement is very good. This result shows the 
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importance of evaluating €j to determine whether the theory can reasonably be applied 
to any particular configuration and free-stream conditions. 

As an additional comparison with the present method and the experimental data 
drag coefficients based on the simple Newtonian pressure equation Cp = 2 sin^^ and 
on inviscid conical flow were calculated and are presented in figure 5. The Newtonian 
prediction and inviscid conical solution drag values are low since viscous interaction 
effects on the surface pressure become important on high-fineness-ratio bodies at high 
Mach numbers. 

As noted in reference 22, entropy layer effects become important for power-law 

y + i 

exponents less than m (m » 0.63 for y = 1.4); therefore, the theoretical predic- 

2y + 1 

tions (which do not include these effects) can be expected to be poorer in that range. A 
less subtle limitation occurs at m = 0.5, where = 0; that is, the ratio of shock-wave 
radius to body radius becomes infinite. This case is the "blast wave" solution for blunt- 
nosed bodies of negligible thickness (for example, a cylindrical rod) as described in 
references 15 and 22. For bodies with nonzero radius, as in figure 5, the predicted 
shock-wave radius goes to infinity as m - 0.5 and so does the wave drag. Thus, the 
theory is not useful for the blunter shapes. 

Wing, conical-body configurations .- Theoretical estimates fpr wing conical body 
configurations can be compared with the experimental data in reference 23. The bodies 
in this reference were halves of right, circular cones, corresponding to m = 1. The 
wings were thin flat plates. The normal- and axial-force coefficients for configurations 
with the first-order Mach number and boundary- layer thickness corrections to the wing 
planform shapes are presented in figure 6. The present theory is in good agreement 
with the experimental data near an angle of attack of 0°, but deviates from it elsewhere. 

The deficiency in the angle -of -attack method is such that the errors in C A and C N 
are generally about equal and in the same direction. This condition results in the good 
prediction of the lift-drag curve (drag polar) shown in figure 7, which produces lift -drag 
ratios agreeing well with the experimental values. Figure 7 also shows the pitching- 
moment coefficient, the theory generally agreeing well with experiment near a = 0°. 

Other data for comparison with theory may be found in references 6 and 7. The 
wings for the configurations tested had delta planforms with several leading-edge sweeps. 
Consequently, they cannot match the shapes used by the theory, but at small angles of 
attack, where the wing alone produces little lift or drag, the aerodynamic coefficients 
should be comparable if they are based on the areas of delta wings approximating the 
theoretical planforms. Figure 8 shows such a comparison at Mach numbers 6.86 (ref. 6) 
and 20 (ref. 7). The theoretical drag polars and the lift variation with angle of attack 
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at M = 6.86 agree well with the experimental points, especially near a = 0° (fig. 8(a)). 
The pitching moment about x = 2i/3 is predicted well near = 0, but the slope shows 
an almost neutrally stable trend whereas the experimental data show the configuration to 
be somewhat more stable. The difference in the distribution of wing area between the 
experiment and theory would contribute to this effect. 

Maximum lift-drag ratios for the same configurations (and some with smaller cone 
angles) in helium at Mach 20 are shown in figure 8(b). The predicted values agree fairly 
well with experiment considering the differences in wing shape and area. 

Example Application of Method 

The preceding comparisons with experimental results have shown that the present 
theoretical method gives good predictions of the lift, drag, and lift-drag ratio and fair 
estimates of the pitching moment for small angles of attack as long as the basic assump- 
tions of the theory are met. Thus, the method should be useful for studying the general 
characteristics of the power- law-body flat-wing configurations at high Mach numbers. 

Just two parameters, the power-law exponent m and the fineness parameter f, com- 
pletely specify these body shapes. For the wings the Mach number is the principal addi- 
tional parameter required, although the Reynolds number, ratio of specific heats, and wall 
temperature also enter through the boundary- layer growth perturbation. In order to 
assess the effects of these three main variables, the theory was used to predict the aero- 
dynamic characteristics of a family of full-scale configurations at two Mach numbers. 

The chosen altitude was 30 km for which the unit Reynolds numbers are 2.21 x lO^/meter 
and 4.42 x 10®/ meter at the chosen Mach numbers of 6 and 12, respectively, based on the 
1962 standard atmosphere (ref. 24). The body volume was set at 2500 m^, giving lengths 
of 28.2 m to 78.1 m (approximately 92.5 ft to 256 ft), for 0.63 = m = 1 and 2,5 = f = 10.0. 
Additional assumptions were y - 1.4 and a ratio of wall temperature to total temperature 
of 0.41667. For each Mach number the range of the fineness parameter was chosen to 
keep 6^ « 1 and « 1. The results of these calculations are presented in figures 9 
and 10. 

Effect of power-law exponent .- Varying the body power- law exponent while holding 
the fineness parameter constant at f = 5 for Mach 12 flight at an altitude of 30 km pro- 
duced the curves shown as figure 9(a). The drag polars in the range 0.63 = m = 0.75 all 
cluster together, and hence so do the lift-drag ratios. Only in the conical case (m = 1) 
does the drag fall significantly higher and the lift- drag ratio lower. The pitching moment 
does show a major variation with m, both in slope and intercept. As m decreases 
from 1.0 to 0.63, that is, as the nose becomes blunter and the aft end less flared, the 
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zero lift pitching moment C mj o increases. At the same time, the stability decreases, 
configurations with m < 0.75 becoming unstable for the moment reference center at 
x = 0.6f,y = 0. 15?^ g. Since the effect of m on the lift -drag ratio is relatively small, 
this parameter could be chosen to minimize the trim drag. 

Effect of fineness ratio .- The computed characteristics for a range of values of the 
fineness parameter f are shown in figure 9(b), For this family of configurations, the 
power-law exponent was set at m = 0,75, and the curves are for Mach 12 flight at 30 km 
as before. At low values of f the peaks in L/D are low and broad and become higher 
and sharper as the bodies become finer. The stability of the configurations is practically 
unaffected by variations in the fineness parameter, as indicated by the almost parallel 
pitching- moment curves. 

Effect of Mach number .- Figure 9(c) shows a comparison of Mach 6 calculations 
with those for Mach 12 for configurations having three of the power-law body shapes. 

Note that the change in Mach number makes a change in the wing planform for each body 
shape. The effect on the drag polars shows clearly in the three sets of curves. At 
Mach 6 the zero lift -drag coefficient Cp^o is higher but the drag due to lift is lower 
than at Mach 12. Since the curves cross before (L/D) max is reached, the Mach 6 
curves of L/D peak higher and at larger Cl values than the Mach 12 curves. If 
the same reference area had been used, the Cl difference would have been larger 
since the Mach 12 design wing is smaller. The pitching -moment curves are little 
affected by the Mach number change. 

Summary of calculations .- The results of the Mach 6 and 12 calculations for flight 
at 30 km are summarized in figure 10. As was indicated in figure 9, the effect of the 
power- law exponent m on (L/D) max is relatively small. For the low fineness ratios, 
the curves form broad maxima centered near m = 0.7; they become more peaked and 
move toward m = 0,8 as the fineness ratio increases. This result compares with the 
value m = 0.75 determined from the Newtonian pressure law as the power-law exponent 
for minimum drag bodies under length and diameter (that is, fineness ratio) constraints. 
There is a stronger dependence of the associated lift coefficient C L (l/d) on the 
value of m, particularly for the less fine bodies. The effect of the fineness parameter 
on (L/D) max and ( -'L,(L/D) max is opposing in that increasing f increases (L/D) max 
(and its dependence on m) but decreases C Lj (L/D) max (and its dependence on m). 

(At any given lift coefficient in the range of calculation, however, L/D can be increased 
by going to a finer body; see fig. 9(b).) The curves of «(l/D) are included in fig- 
ure 10 in order to show that the calculations of (L/D) max occur within the range of small 
angles of attack for which the present method gives its best results. (See fig. 6.) 
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CONCLUDING REMARKS 


This paper has presented a method for calculating the longitudinal aerodynamic 
characteristics of a family of configurations in hypersonic flow. These configurations 
each consist of a half-axisymmetric power-law body surmounted by a thin flat wing for 
which the planform matches the analytical shock-wave shape about the body at an angle 
of attack of 0°. The method is based on the power-law similarity solutions of the hyper- 
sonic small-disturbance equations. These solutions require three basic assumptions: 
the Mach number is large, the body is slender, and the shock wave is strong. A first- 
order perturbation allows the calculation of Mach number effects, and a perturbation to 
the body shape provides for the boundary- layer growth. Skin friction is accounted for 
by using compressible, laminar boundary -layer solutions at the computed pressure dis- 
tributions integrated over the body and wing surfaces. A computer program has been 
written implementing this method; sample computations using the program have taken 
only a few seconds per case. 

When compared with experimental data for axisymmetric power -law bodies and 
for wing — conical -body configurations, the present method gave good agreement where 
the basic assumptions were satisfied. An example series of computations with varia- 
tions in the principal parameters at a full-scale flight condition showed that varying the 
power- law exponent has a greater effect on longitudinal stability and trim than on the 
lift-drag ratio. The computations for Mach 6 gave higher maximum lift-drag ratios, 
higher drag coefficients at zero lift, but essentially the same stability characteristics 
as their counterparts for Mach 12. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 25, 1973. 
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APPENDIX 


COMPUTER PROGRAM FOR CALCULATING 
THE AERODYNAMIC CHARACTERISTICS OF POWER- LAW 
WING- BODY CONFIGURATIONS 

The calculation procedure described in the main body of the paper for obtaining 
the aerodynamic coefficients for power-law wing-body configurations at hypersonic speeds 
has been programed for high-speed digital computation. The program will also compute 
the zero angle -of -attack drag for an axisymmetric power-law body alone. The purpose 
of this appendix is to provide a description of the necessary input and available output as 
well as a FORTRAN IV (ref. 25) listing of the source program. A separate program to 
compute the parameters needed to obtain two input values from the figures of reference 20 
is also listed and described. 


Description of Program 

First, the program reads all the input variables describing the case to be computed. 
After calculating geometric constants, it goes through the angles of attack, computing the 
body axis forces and moments, interpolating the similarity solution parameters from a 
stored table. Skin friction is calculated for each angle of attack and added to the axial 
force. The results are then transformed to the stability axes. If at least three angles 
of attack are included in a case, a quadratic interpolation of the drag polar is made to 
obtain (L/D) max and other quantities, which are printed out along with the body- and 
stability-axis coefficients. A summary subroutine assembles certain quantities for sep- 
arate printout after completion of all cases. 


Program Listing 

The FORTRAM IV listing of the source program used on the Control Data series 6600 
computer system at the Langley Research Center is as follows: 


PPCGPAM HYPAERCM INPUT *201 t OUTPUT =401 1 T APF 1 '® INPUT * TAPE 7*601 ) A 1 

C 

C HYPERSONIC AERODYNAMIC CHARACTERISTICS OF POWEP-LAW WING-BODY CONFIGURATIONS 
C 

DIMENSION HE AD I 8 ) » Y(6>i YE(*5)» VAPD(I3,fc), VAPI(13)» ANGL(ll)t A 2 

1 ST NE ( 1 1 ) t COSE Cl J If PBBP0UU), CFDCFOldllt PUPIOIllDt CNBdllt A 3 

2 CNDCllJt CNWO.ll, DP5RY0S dll* CNdl), CAP(U), CAFClll t CAdl), A A 

’ C.KIMt COdl I, CLCOClUt CMCGU1), CDANGL(11.2lr CDAL0I2I, A 5 

* CDALNI 2) » X ( 19 1 » XW(19), PB(19I, DSWCCSU9), DSWCSL(7,*>, A A 

s rcssinei, txse(I8>» delr<18I, ptwpui8> a 7 

EQUIVALENCE (Y,ETAB), (Y(2I,F0), <Y(6),A1), «YE,ETA6E», a 8 

1 (YEI2),F0E>, (YE(3),F1E), ( YE (4 I f 0 JOE I ♦ I YE I K I »0 JIE » , A 9 

P I ANGL »CO( 12 1 1 COANGLI 12)1 A 10 
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APPENDIX - Continued 


COMMON &P,Z1ME,E2U2, FOE, F1EK»THG12G,TWM,EMBF2,02ME2K,X2M, WING, F01, A 11 

1 Fhk,DL9»ETB»CNETB,XSE,GMA,GM1,GM12*GP12,EM,EMT,EM32>ZMM,ZEMI1» A 12 

2 THP2M, AMCH2,0EL,TWEM,AEP A 13 

NAMELIST /OATA/NCASE,GAM,TlNF,AMCH,NALF, ANGL » EM, F ,PEL ,SSB, XCG, YCG, A 14 

1 PP.P T * C6 FACTR* ANGO, XOUT , PB 8P0L »CFDCFO l , ALAMf R , BDYONLY*DEL, EPStXSE A 15 

2 /CUT/T,ANG, ANG0,FME,YE,EMI,CAFR2,CAFU,CAFL,AI,PTWP1 A 16 

r VART = fm, VARO = ASSOCIATED VALUES OF ETAB,F0,F1,DJ0»DJ1,A1 FROM TABLE I 

DATA X/O.,. 0003,. 0006,. 0009,. 0012,. OOIS,. 0024,. 0036,. 005,. 0085, A 17 

1 .CIS, .025, .045, .08, .14, .25, .45, .7,1./, Pi ,P!F,0T0R/3. 14159265359, A 18 

2 2.6586808, .01 7 45329252/, KP ,L LI M, ML , ML1, MCA SE/-1 , 19, 18 , 17, 1 / , A 19 

3 VART , VAR 0/1 .,.95, .9, .85, .8, .75^.7, .666667, .633333, .6, .55, .53, .51, A 20 

4 .914034, .91034, .90465, . 89743, .88798, .87507, .85648, .8388,. 81391, A 21 

K .77647, .66414, .56901, .37221, . 87445 ,. 84711 ,. 81 63, . 78174 ,. 74265, A 22 

6 .69806, .64662, .4^763,. 56403, .51478, .42678, .385, .33757, .9179, A 23 

7 1 <0591,1. 23 06 , 1.4?86, 1.6897,1.9811,2. 29 86 , 2.4964,2 . 6392 ,2.6593, A 24 

8 2 .2*! , t .8876,1.411, .07323, .07589, .07909, .08303, .08799, .09444, A 25 

9 .10318,. 11098,. 12129, .13564, .17318, .20119, .25443, . 0841 ,. 09551 , A 26 

X .11044, .13067, .15918, .20098, .26458, .32562, .40794, .51763, .74598, A 27 

1 .86687,1.0411, . 47549, .52709,. 58604, .6 e 291,. 72741,. 80732, .88631, A 28 

2 .932^6, .96566, .98034,. 96791, .96377, .97539/, GAM/1.4/, A 29 

3 SSB,TWTT,PBP!,CAFACTR, XE,xnuT/0.,. 416667,1. ,1.,. 10, .FALSE./ A 30 

LOGTCAL TFMTN,8DVONLY,WTNG,XOUT A 31 

EXTERNAL FUN1 A 32 

RDYONLY*. FALSE. A 33 

1 °EAD 38, HEAD A 3 * 

IF ( ENDFI L c 5) 32,2 A 35 

2 RF AD DATA A 36 

CL DMX=0 . A 37 

COMT N=1, . A 38 

1FMI N=. FALSE. A 39 

GO TO 12,4,5,6) , NCA.SF A 40 

C, NCASE = 1) GAM, 2) AMCH.TINF, 3) EM, 4) F,PEL,SSB HAVE NEW VALUES 

3 GM1* GAM-1 . A 41 

ZGM1=2./GM1 A 42 

GM12*. 5*GM1 A 43 

CPl*G*M*l. A 44 

GP!?=.5*GP1 A 45 

F01=1./GP12 A 46 

GM A = G.AM A 47 

GPJ4=.28*GP1 A 48 

ZGG1=GAM*ZGM1 A 49 

ZG!G=-4./ZGG1 A. 50 

GM 1=1. /GAM A 51 

ZG2G=GMT-. e A 52 

GP1G=GP’*GMT 6 53 

GF*.12*<?.*GAM)*#1 .S*GP12**GP1G A 54 

GX1***GM1*GMT a 55 

G’7=-. ">7/6X1 A 56 

THG1?G=1.-.5*GX1 A 57 

RT8GT»1./S0RT(8.*GAM) A 58 

AU = 1.2!7*TwTT+.4704 A 59 

AT2-TWTT+.246 A 60 

RTTWTT=S0RT<TWTT) A 61 

4 AMCH2 = AMF H*AMCH A 62 

AMCH3= AMCH*AMCH2 A 63 

4MI=1./AMCH2 A 64 

TT8TI-1.+GM1?*AMCH? A 65 

pTTT8TT = SQRT(TTBTI ) A 66. 

SUTH = !98.6/<TINF*TTRT1 } A 67 

4LAM= ( (l.+SUTH) /( TWTT+SUTH I ) *RTTWTT A 68 

Sijtht = 1 96.6/TTNF A 69 

C2*( C 1.4-SUTH) / ( 1 . + SUTHI > )*RTTTBTI A 70 

PTC1=SQPT ( AL4M/C2. J A 71 

tpTT=.273+(.19 c +.53?*TWTT)RTTBTT A 72 

pTC3 = S0PT ( ( l.T’SUTHl \( ( TPTT+SUTHI I* SORT f TPTI 1 ) A 73 

GCMG=GE*PTC1*AMCH**GP1G A 74 

P0=2.*AMI*GMT A 75 
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APPENDIX - Continued 


5 FMPl -FM+l # A 76 

*=MP 123.f5*PMPl A 7? 

f=MP2 = FM + 2. A 78 

F^H1=CM-l. A 79 

7MM1=?.*FMM1 a BO 

TWP=?.*FM A 81 

THP2P=8.-TWM A 82 

ZMN=?.-FM a 83 

EMI=1./FM A 84 

ZFMIl=?.*<EMT-l. ) A 88 

FM32=l.5*EMT-2. A 86 

ZGH!32=i .5*ZNM1*GMI A 87 

ZMP1I-1 ./(EMP1 + EM) A 88 

TWN3I=!./(FMPl4FMP2| A 89 

FM*I=t./(4.*FM-l.) A 90 

SIXMI=i ,/(6.*EM-l. I A 91 

THRMMI=1./I3.-EM| A 92 

FIV2M!*1 ./C'.-TWM> a 93 

FM4019=.4019*EM A 94 

FMX1!*] ./(1.8303-FM4019J A 95 

FMX152M=FIV2MVEMX1I A 96 

EMX2=.5/< .9274+EM4019) A 97 

TP=20.*<1 .Of— EM) A 98 

AI=AT1-AI2*£MM1/(EMM1-G37/EM4I ) A 99 

CALL MTLUP « EM • Y , 2 ♦ 13, 1 3 , 6 , T P, VARI , VARO I A 100 

ET ABI = 1 »/ETA8 A 101 

0NETET=ETABI-1. A 102 

ZG1ET=ETABI/GP12 A 103 

ETp=ETAB A 104 

0NFT8=1./C1.-ETABJ A 105 

AMCHN=(l.-ZMMl*CAMCH-l.n**GMl A 106 

GETIEMF=GCMG*ETAB*AI*SQRT(EM4T >*(EM*EM*F0»**ZG2G/AMCHN A 107 

F11=F01*( A14(ZGGl+2./GPlG-FMT/GHl)-EPI*EM!/ZGGl) a 108 

6 FT*!./F A 109 

EMBF?=(EM*FI)**2 A 110 

FSO=F*F A 111 

FM1M2F2=ENP1*EMP2/FSQ A 112 

PTPELI=J«/SQPT(REL) a 113 

ES6V=GM1*AMCH*FSQ*PTC1*9TRELI*RT8GT a 114 

0EL=ET4BT*F! A 11* 

FPS=AM»/10EL*DEL» A 116 

A«=P-41*EPS A 117 

EMSA V=4.*ESAV*AT*ETAB/(EM*SQPT (F0/EM4 I > I A 118 

TWEM = EMSAV*.S ■ a 119 

PSE=SQnT(GETTEMf*PTPELl*DSL**ZGlG*XE**ZGM132> A 120 

PT2P1=PT0T(PSE) A 121 

RTPOLt=SORT< AWCH*PTTT8TI**3/(C2*PEL*PT2P1)) A 122 

AFPX*AFP*EMX1F2M A 123 

5M22X=.25*EMSAV*EMX1I A 124 

G?PRr=.69053*GM12*AMCH3*RTC3*RTPELI A 125 

EMFSAV=0. A 126 

WING GEOMETPir. PARAMETEPS ANO FLOW CONSTANTS 

SISTSB=ETABI*( 1 . *FMP1* ( AEP*THP WM1+EM S 4 V*F 1 V2MI) I A 12T 

SBBYS=1./SISTSB A 128 

IF (SSB.GT.O.) SBBYS=1./SSB A 129 

S7ST8YS=S1ST$B*SBBYS A 130 

OEL?=DEL*OEL A 131 

PtMS=.5*PT*EMPl*SBBYS A 132 

ot FPFL! = SQRT(C1,+AEP«-TWEP)/(P1F*(FMX2*AEPX*EM32X)I> K 133 

XSC=ENP12 a 134 

ALAMCR=G2MRC*PTCREL! A 135 

POL AW=-po*ALAMCR A 136 

CF0CR=1.328*RTC1*PTPELI*RTCRELI*S1STBYS A 137 

CAFLC-2.*ALAM*AMI*RTR0LI A 1*8 

CFCCN=PT2P1*C4FLC A 139 
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APPENDIX - Continued 


FMET$=EMPl*ETABt*$BBYS A 1 *° 

M =3 A 141 

XK 1 =(ML- 1 )/M +1 a lA2 

kk=kkj.+i a lA3 

K=KK A 144 

MK=M*U-2)+l A lA5 

OSWroSdLlM} = PTM$ A 1AA 

PB<LLTM)=ETAB a !47 

LKl I» A 14 ® 

C 0 SSL = 0 . A 1A<? 

PS= 1 .+TWEM+AEP A 150 

XW(KK)=| . A 151 

C 0 SSUKK)= 1 ./SQ»T< 1 .+EMBF 2 > A 152 

00 8 11=1**11.1 A l 53 

L=LUM-LL A 1*4 

Xt = X ( 1 ) A 1 5 5 

XM = Xl.**FM A 156 

0 SWC 0 S(L)=PIMS*XM A 157 

IF ( FSOYONLY) GO TD 8 A 158 

PR (t )=FTAB*XM A 159 

IF (MK.NE..U GO TO 8 A 160 

K! =K A 161 

K=K- 1 _ 4 162 

MK=M* ( K- 2 ) +1 A 163 

XW (K I =XL A 164 

XSF=. 5 *(XW<K 1 »+XL> A 165 

TXSEIK)=XSF A 166 

PT.taPl C K 1 =PT 0 T 2 (RSEI A 167 

DFLP{K>=PSE-ETA 8 *XSE**FM A 168 

RSK 1 =PS A 169 

XM 1 ? = (XM/XU **2 4 170 

PS=XM*( 1 .+<TWE!VSGPT<XL)*S 5 P)/XM 12 ) a 171 

C 0 ?SUK|« 1 ./SQRT( 1 .*EMBF 2 *XM 12 ) A 172 

0 RM 1 S=EMETS*CRSK 1 -RS-R 8 (L 1 »+RBCU ) A 173 

DShCSUK t K)=DRMl$*. 25 *(raSSL(K)+COSSL(Kin A 174 

Ll=t A 175 

CO 7 J*K 1 ,KK A 176 

OSWOSU J,K> = DRM 1 S*C 0 S 5 LUI A 177 

CONTINUE A 178 

05 WC 0 S=PTMS*(. 1 *X( 2 »)**EP A 179 

I F (RDY 0 NLY 1 GO TO 10 A 180 

XW= 0 . . A 181 

XSF=.R*XW< 2 > A 182 

TXS£=XSE A 1®3 

PTWP 1 =PT 0 T 2 (PSE» A 184 

OELP=PSF-ETAB*XSE**EH A 185 

0 RM 1 S=EMETS*(RS-RB(U> > A 186 

OShrSL=DR»US*. 25 *COSSl ( 2 > A 187 

DO 9 J= 2 »KK A 1 B 8 

DShCSI < J ,1 I=DPM 1 S*C 0 SSU J) A 189 

XLl = xm A 190 

XL 2 =XC 16 ) A 191 

Pll=R 8 ( 6 > A 192 

P 12 =R B ( 1(1 A 193 

ANGIE OF ATTACK VARIATION 

#NGO=OEt*(RB(MLI-FTAB>/tl.-X(ML)> A 194 

10 PRINT DATA A 195 

PRINT 36 , HEAD A 196 

NANG=NALF A 197 

DO 28 T = ! , NANG A 198 

ANG = DTOP*.ANGl ( I ) A 199 

IF (NOA 5 E- 3 ) 12 , 11,11 A 200 


19 



OS 


292 v 
192 V 
092 * 
tn * 
862 * 
162 V 
962 V 
662 V 
*62 V 
£62 V 

zn * 


162 * 
062 V 
bn v 

8*2 V 
1*2 V 
9*2 V 
6*2 V 
**2 V 
£*2 V 
2*2 V 
1*2 V 
0*2 V 
6£2 ? 
8£2 V 
1£2 * 
9£2 V 
S £2 V 
*£2 V 
££2 V 
2£2 V 
t£2 V 
0£2 9 
622 V 
822 V 
122 V 
922 V 
622 V 
*22 V 


£22 V 
222 * 
122 V 
022 * 
612 V 
812 * 
112 V 
912 V 
612 V 
*12 V 
£12 V 
212 V 
112 V 
012 V 
602 V 
802 V 
102 V 
902 V 
602 V 
*02 V 
£02 * 
202 V 
102 V 


»a*u=j-»nj 

Mdfcd id- >13 I d 
AV0*Sdd=Xd 
w l>/d* t • 1-fcdWZ ) /* 1=A *S3d 
i ldhz*2n*d3=wn*d 

( m3W dd*3Od)ia0$*23W3)/3W2dHi* 38 tfid^IVsAVS 3 U2flV 7=207:13 

S^20“= JS«2Q 
SW2G** Z=2SrlcQ 
SAbd$*>tcdrtZa = SrtZa 
A* J/23rt3*21d0 = >l2342Cj 
2**ltf di=j»2n*AV J 

S38nSS3ttd AG09 CJNV ONIH G31Vd031Nl HOdd SlNdijIdddGO SiXV AGGb 

iJidNG*ii3id^3oru=bidra 
30d*l drt3=lWdGd 
Wd*3Gd=Kd0i 
23*d/i*0 = l2d*‘J 
3H7-*2=dWtZ 
( 13h bd*l£3-I 3w3 )/ldW3*2i^-ll9-i? 

•I-3*3=X3w3 
3rt2-*£-8rt2dHi 
*d+dwZ=HdftZ 
uWd* * 2 = 3w2 
* 1 "*3W dd = LdWtd 
dWd*'*=3rfbd 
dW3*dirtb = Z:il*d 
<38Vi3-*tW i = i 1313 
i9Vid*3£7ld=lVdld 
(CldVA 4 iaVA 4 di 4 6*£I 4 £I 4 2 4 iA 4 3W3) QflllW VIVO 

3W3=A^SdWd 
22 *£2*22 (AV$3W3-3*3) di 
0 * L=dWd 
12 01 GO 
16*=dwd 

U 4 U 4 6t Ub # ~dH3) di 
( (Nl*nd-nx ) / tN 1*21 d-21X ) >OOT9/<lVdOhd)OirW=dhd 

8i 4 61 4 6t UVauHd) di 
( Nl* nx+11dj / fNl*21X*21d) = iVciOHU 
UZ G1 00 (l-Od^d) di 
Z**lNl*d+*X)/* l-m 
( 1 13S00/U >*im»Ni. 

XjViiV dD 310NV 01 300 dcVHS AUD8 3Alij3dd3 

(liliOIdfld-l in0d98d^vfV10dUSA8iSiS«( iJSOAbSdG 
< lUDldOd 4 i DION* 4 *Z lfslldd 
i -I-dOd)*Od=< i) 101 dOd 
*G=30d ( # 0 B lV^o) dl 
I0D7**<0*3*2Iw0-*I)=3Gd 
12 01 DO 
Wd=3Wd 
1=01 
* l=zn 

U)10d96d*WVlGd*SA8lSlS=< 1 JSOAflSaU 

be Oi UO 
( I nONV 4 6t INldd 
1-0N9N=0^VN 
11 Oi 00 (ONVdVOflNV) dl 
(11 lOIdlldH D'ttNV 4 ££ iNlbd 
t * t-SOd ) *Od=l 1 > lOIdOd 
( (2***0>4ldO+ # I ) 1 dOS+*UXl dO ) *0^ VO* * I = S Gd 

CiV3«*4dO = 70XldO 
9l 4 ^i 4 £l (0 NV> di 
ONV*HOWV)$av-JVO 
(2**< i )3NlS-*Uidt>6 = l DdSOO 
<0N* IN 1 S “ I i ) dMIS 
li 4 6l 4 *l (ON*) dl 
ddHSSddd HOVddflS odddfl OiUrt 


12 

0 

J 


22 

12 

02 

61 

8 i 


U 

J 

11 


9 1 


61 


*1 


£1 
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APPENDIX - Continued 


t 

r 


r 

c 


24 

C 

r 

25 


26 


27 


28 


C 

r 

c 

29 


EKGO^FK*GME' , T 
FP$AV0=F1FK-P8PI*FKG0 
FPSAVI=F! FK-EKGO 

FP5AV2=CINFTFT*( FK*0 J1F*FTE1T— FKGO ) 

CNP(T) = D2MS2*(F0EM1*FFSAV‘*EPSAV1 ) 

CND ( I>=D2MS?*FOFMl*EME*ETP AT*EF4U2/< < ZME + l • I *THR 2ME ) 

CNVM I )=D2MS2+< FHPl*t DJETB*EFSAV4*ZGlETMEMSAV/( FRME1-ZMM1 l«*AEP/ 

1 i ) )-EMF*ETPAT*E4UM*F0F/THP 2PE)+F d SAV 2 ) 
f AP{ I )=PT*D2MS*EMP1 *F f *( FOEM*<1. . / 1 ZME +ZMH1 >+FE4U2*FM4I )+.5*EPSAVO) 
CNt I )=f NBC I ) ♦TNDC H +CNW< T ) ♦DPSflYOSM > 

rMN*02MSC* p MPl*( FMP2*( (FOE+DJETB ) *< l./ZMEMff E4U2*TWM3I ) + 

1 ZGlET*(*MSAV/( FPMF+THR2M) +AEP/< ZMEM-ZMM1 ) M +EPS AV1+FPSAV2 ) 
CMA-D2MSC*EM1M2F2*< FOE** (1 ./( ZMFM + ZMMl ) +FF4U2*S IXMI J+EPSAVO/3. > 

r^rr, ( t )=cmn+cma-xsc*dp$pyqs+.5*fmp2* <cnu ) *xcg+c ap< t )*fi*ycg> 


SKIN FRICTION ON BODY 
FOF^FCF 
E2U2=.5*E E4U 2 
&PC=GAM*F*E 2 /FK 
«P=APO/PT2Pl 
WTNG=. FALSE. 

CAF82 = CFCON*SKNFPC ( X, LLI W f OSWCOS,GAM , *L AM, FUN1 I 
IF ( BDYONLY ) GO TO 2* 


SKIN FRICTION ON WING UPPER AND LOWER SURFACES 
CAFU^CFOCR* ( 1 . + CFDC FOl ( I ) ) 

C AFL=0 . 


Wl NG= .TRUE . 

DO 24 K-I,KK1 

Kl-KK-K+1 

AP«APO/PTWPl(K) 

DtR-DELR ( K ) 

S A VXW=XW ( K ) 

XW(K)=TXSF(K) 

C AFL=CAFL*PTWP1<KJ*SKNFRCW(XW( K) , Kl, DSWCSL(K t K) f GAM, ALAM f FUNl > 


XW | K ) =SA VXW 
CAFl=CAFLC*CAFt 
IF ( XOUT ) PRINT OUT 


LIFTf OPAGt ANC L / D 
CAFU)=CAFB2+CAFU+CAFt 
GO TO 26 
C AFI T)=?.*CAFB? 

CAP(I1*?.*CAPUI 
CN-0 • 

CMCG-0. 

C A (T )=CAP( I ) +C AFACTP*C AF (I ) 
cun-c Nm*cosEtn-cAm*siNF(i) 

COm-CAf I l*COSFm+CNU l*SINE(I) 

neon >=clu wco< n 

IF (CDMIN.LT *CD Ml) GO TC 27 
I MIN- 1 
r DMIN-COU ) 

IF (CLCDU KLT.CLDMX) GO TO 28 
CLCMX=CLCOn ) 

TMAX=I 

CONTINUE 

IF (XOUTJ PRINT 36, HE AO 
IF ( NANG-3 > 31,29,29 

QUADRATIC INTERPOLATION OF DRAG POLAR to get (L/D)MAX, etc * 

IF ( JMAX.LT. 2) I MAX=2 

IF (IMAX.GE.NANG) TMAX=NANG-l 

IXP=IMAX+1 

T XM-TMAX— 1 

yi=cd<ixm) 

Y?-CO(TMAX) 


A 263 
A 264 
A 265 
A 266 
A 267 
A 268 
A 269 
A 270 
A 271 
A 272 
A 273 
A 274 
A 275 
A 276 


A 277 
A 278. 
A 279 
A 280 
A 281 
A 282 
A 283 


A 284 
A 285 
A 286 
A 287 
A 288 
A 289 
A 290 
A 291 
A 292 
A 293 
A 294 
A 295 
A 296 


A 297 
A 298 
A 299 
A 300 
A 301 
A 302 
A 303 
A 304 
^ 305 
A 306 
A 307 
A 308 
A 309 
A 310 
A 311 
A 3X2 
A 313 
A 314 
A 315 


A 316 
A 317 
A 318 
A 319 
A 320 
A 321 
A 322 
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Y?=CD(TXP) 

xt=rutxMi 

X? = CU IMAX) 

X?=CL ( I XP ) 

X12=X1-X2 

X2?=X2-X? 

TF ( I FM T N ) GO TO ?0 
X 3 1 = X3— XI 

A=<Yl*X23+Y2*X3HY3*X12)/(-X12*X23*X31> 

XA*.5*< A*(X3+Xl)-< Y3-Y1 >/X31 ) 

XA2YAMY2 + X2M 2.*XA-A*X2 ) 

CLMX2=X A2YAA/A 

IF { CIMX?.L T * 0 * I GO TO 30 

tLMX=SQRT(CLMX2> 

CDMX=2.*(XA2YA A-XA*CLMX) 

CLCMX=CLMX/CDMX 

TAIL MUUP <0.,CDAL0,2,NANG t lX f 2,KP,CL,CDANGU 
TALL MIL UP (CLMX,ALPHX,2,NANG,ll,1 , I M AX , CL , ANGL > 

TFMIN^.TPUE. 

IMAX=1MIN 
GO TO 2° 

30 Y322!-(Y3~Y2) / ( Y 2-Y1 ) 

X3221=X23/X12 

CLHN = .^*(Y322TM X2*Xl)-X3221*< X3 + X2) ) / i Y32 21-X3221 ) 

CALL MTLUP CCLM^ f CDALN,2,NANG f 11, 2 ,KP , CL ,C DANGU 

r 

C MATN OUTPUTS 
C 

PRINT ?7 f CLOMX, ALPHX,CLMX,CDMXtCDALC,CDALN,CLMN 

CALL SUMMARY fCLDMX f ALPHX, CLMX # COAL 0 , C DALN, CL ( 10 ) « HEA C, NCASE ) 

31 PRINT 34# (ANGUI)tCLU)tCOm f CMCGUI»CLCO(T) f CNm,CNBtn f CNO<II 
IfCNWn J,nPSRYOS(T) t CAU)fCAP(n f CAFU>#T*ltNANGI 

NCASE=4 
GO TO 1 

32 CALL PPNTSUM ( CLDMX, ALPHX, CLMX, C DA LO , COAL N,CL« 10) , HEAD, NCASE ) 

STOP 

C 

33 FOP MAT <F12.2,15H DEG* PUPIQI ^,F10.5) 

FORMAT <//3X,5HAlPHA,8X,?HCL, AX,2HCD»8X,2HtM,7X f 3HL/D,9X,2HCN,8X f 

1 ^FCNB,7X,3HCNO,7X,?HCNW f 6X,4HOP/Q,8X, 2HCA, 8 X, 3FC AP , 7X,3 HCAF // 

2 (F7 .1,3X,3F1C.5,F9.2, 2X,5F10.5 f X,3F10.5n 

35 FORMAT <//F 8 . 2 ,?E!H DEG., TOO NEGATIVE FOR BOOY) 

36 FORMAT < 1H1/ 20X 8A1 0 / ) 

37 FORMAT C//T 1 H <L/D)MAX =,F8„4,11H AT ALPHA =,F7.4, 26H DEGREES, 

1 with cl and CD =,2F9.6//6H CDO =,F1C.8,14H, AT ALPHA 0 *,F8.4// 

2 9H CD MIN =,F10.8,l?H f AT ALPHA - t F8.A,9H AND CL =,F 8 . 6 ) 

38 FORMAT ( 8A10 ) 

FNO 


A 323 
A 324 
A 325 
A 326 
A 327 
A 326 
A 329 
A 330 
A 331 
A 332 
A 333 
A 334 
A 335 
A 336 
A 337 
A 338 
A 339 
A 340 
A 341 
A 342 
A 343 
A 344 
A 345 
A 346 


A 347 
A 348 
A 349 
A 350 
A 351 
A 352 
A 353 
A 354 

A 355 
A 356 
A 357 
A 358 
A 359 
A 360 
A 361 
A 362 
A 363 
A 364 
A 365- 


C 

C 

C 


c 


1 


FUNCTION PTOT <RS) B 1 

TOTAL TO STATIC PRESSURE RATIO ACROSS SHOCK, AND SHOCK POSITION B ^ 

COMMON DUM9<16], XS , GAM , GM1 , GM1 2 , GP12 , EM , EMI , EM32, ZMM, ZEMU , THR2M, B 3 

1 AMCH?,DEL # TWEM,AEP B 4 

X POSITION FOR GIVEN SHOCK RADIUS B 5 

XSMO-PS B 6 

X0=RS**EM! B 7 

DO 1 T = 1 , 10 B 8 

X01M=X0**EMI B 9 

XSM=P S/(l.'M4£P+TWfH/SQRT(X01M))*(X01M/X0>**2) B 10 

IF ( XSM/XSM0.GT*.999) GO TO 2 B 11 

X0-XSM**EMt e 12 

XSMO-XSM B 13 
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? XS^XSM**EMT 

GO TO ? 

FNTp Y PT 0 T 2 

C SHOCK RADIUS AT GIVEN X POSITION 
X$M=X$**EM 

RS*XSM+( AEP+TWEM/SQPKXS) )*XS*XS/X$M 
3 X$L?M=XSM**2 

TTHX?x(OEl*< EM*X$L 2M/XS+XS *ZMM*A EP+.5*THR2M*T WEM*SQRT ( XS ) ) )**2 
AM?=AMCH2*TTHX2/ ( XSL2M+TTHX2) 

AMG1=GP1 2*AM2* ( 1 ,+GMl 2*AMCH2 ) / < 1,+GM12*AM2 ) 
AMG?*GP\2/(GAF*AM2-GM12> 

PTCT=FXP( ( ALOG< AMG2 ) +G AM* A LOG ( AMG1 ) ) /GHl ) 

RETURN 

END 


8 14 
B 15 
B 16 
R 17 
B 18 
B 19 
B 20 
B 21 
B 22 
B 23 
B 24 
B 25 
B 26 
B 27 


C 

1 


c 


2 


r 


FUNCTION SKNFRC 1 X , UI M, OSWCOS t G AH, At AM, FUNl) 


LAMINAR, COMPRESSIBLE SKIN FRICTION <0ATA ARE FOP TW/TT * ,41667) 


r i l , 


DIMENSION X ( 19 ) » DSWC0S<!9), Xl(19)t U(19)t BETA(19J» 

1 BIT ( 3 ) » BU5). TTH2 < 30 ) » THT2<2) 

EQUIVALENCE { FWPP, THT 2 ( 2 > ) , „ , . _ 

DATA B f TTH2/-.2,-.l»0. ».05,.lt-2,.3».A,.5, .6 t .8 t l.»l.Zfl.6f2.» 

1 , 29 3, . 2A8B, . 2 205 , .209?* . 1 998 , • 1 B45 , . 1725 , .1627, . 1547 , ,1479, . 137 ♦ 

2 .128 c ».!21Bt«!U6,.l'<07, • 269, . 387, .4696 * .5051 5373 t. 59A.A, -^ATt 

— - - * — ---- - * ** f>r\JL /. 6iT7 _ 1 CIA77 .1 _ 1 / . JP/1I. / 


2 ,128 F , -1218, .1116, .1^07, 

3 . 6 °, ,73 165,-7703, .8^08, .9044,-9627, 1 , 0677 ,1-1613/ , 
4 f GX1,GX2,ZGM1/-. 28571*2 8^ 71 43, 1.2857 14285714, 5./ 
EGIN WITH X(l) = 0., XT ( 1 ) - 0-, BETA = .5 (BLUNT- 


BEGIN 

BETA ( 1 ) =• 5 
TMPT2=. 23209 
FPTH=. 28779 


NOSED BODY) 


GO TO 2 

FNTPY SKNFRCW „ 

BEG^N WITH X(l) = XSE, XI ( 1 ) - 0., BETA * 0, (UNOEPSIDE OF WING) 

BETA=G. 

TMBT?=-441 
FPTW=. 22052 
XM-X + -1*< X( 2 )— X ) 


PFPT2“P(XM) 

UM = SQRT(ZGM1*<PEPT2**GXT-1. > ) 

CALL MG A US S < X, XM, T ,XTM, FUNl t FOFX,' ) 

X T M- AL AM*X I M 

TH 2 U 0 =TM 8 T 2 *XIM 

DCF(11=PEPT?**GX2*FPTH*OSWCOS(1I*SQRT(UM**3/TH2UR ) 


XT (1 ) "X I M 


U(1)=UM 

DO 9 I P -2 , LL I M 
XP= X ( IP ) 

CALL MGAUSS ( XM f XP , 1 , OX I , FUN1 , FOFX , 1 ) 

0XI-ALAM*DXI 

XTUP)=DX!+XIM 

pFFT2=P(XP> 

U(IP>=SQRT<ZGM1*(PEPT2**GX1-1.)J 
ALNUP=ALOG(U( IPI/UM) 

ALNUR2=?.*ALNUR 
OXTUTH“DXI /TH2UR 

rLAM“7-7609*(U( IP) /UM— 1 . I/DXIUTH 

88T«PLAM>M1. + (PLAM-1 * >*< (. 020673 7*0X1 UTH-. 20419 ) *OXIUTH+. 344145 ) ) 

KK^l 

K1 =0 

ITERATION FOR BETA (LOCAL VSLOC ITY-VAP I ATI ON PARAMETER) 
no 7 J«1 ,29 
PIT 0 = BET 
PTT(KK)xBET 
TF (KK-3) *,3,3 


C 2 
C 3 
C 4 
C 5 
C 6 
C 7 
C 8 
C 9 
C 10 
C 11 
C 12 
C 13 
C 14 
C 16 
C 15 
r. 17 
C 18 
C 19 
C 20 

r. 21 

C 22 
C 23 
C 24 
C 25 
C 26 
C 27 
C 28 
C 29 
C 30 
C 31 
C 3 2 
C 3 3 
C 34 
C 35 
C 36 

r 37 
C 38 
C 39 
C 40 
C 41 
C 42 

C 43 
C 44 
C 45 
C 46 
C 47 
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3 


4 

s 

6 


7 

8 


9 


10 

C 

11 


KK*1 

P??*BIT(3)-BITf 21 
821=RIT(2)-BTT<1J 

T F MBS< B21+B32 ) .LT.A9S( 832) J GO TO 4 
B0EN0M-B32-B21 
? F [ BnfrNOM • EQ *0 • > GO TO 6 
BFT=(BTT(11*B1T(3I-RIT(2>**2 ) /BDENOM 
K1 "0 
GO TO * 

flPT=<R T Tm+eiTm )/?. 

GO TO * 

KK-KK+Kl 
K 1 =1 

T E (BFT.GT.2. 1 BET*!.* 

CALL F'TIUP <8ET t THT2,2,l5,l*,l f JP,B, TTH?) 

RFT = ALNU^ 2/( ALNUP + ALOGU . +OX!UTH*( 2*-BET) *THT2 ) ) 

TF ( APS ( BET /BITS— 1 . ) — • 0001 ) 8,8,7 
CONTINUE 

PRINT Tit XP , B I TR , RET 
RETA(TP)=RPT 

CALt MTLUP <BET,THT2,2,15,lf , 2, JP,B,TTH2) 

TH2U P = ( 1 •>0XlUTH*{2*“PED *THT2 ) *TH2UR 

DCf ( TP 1=PEPT2**GX2*FWPP*DSWCOS UP) *SQRT<U< TP ) **3*THT2/TH2UR > 

XM^XP 

XIM^XI(TP) 

UM=U(TP) 

TF ( LUM.EQ.2) GO TO 10 

SKNFPC-SUM(X,DCF,LLIM1 

RETURN 

SKKFRC». r *l OfF+OCFI 21 )*(X( 2) -XI 
RETURN 

FORMAT (/26H BETA UNCONVERGED AT X/L =,F8.5,12H BETA VALUES, 2F12*8 
1 ) 

END 


C 48 
C 49 

r so 
C *1 
C *52 
C 53 
C 54 
r 55 
C 56 
t 57 
C 58 
C 59 
C 60 
C 61 
C 62 
C 63 
C 64 
C 65 
C 66 
C 67 
C 68 
C 69 
C 70 
C 71 
C 72 
C 73 
C 74 
C 75 
C 76 
C 77 
C 78 
C 79 
c ao 

C BI- 


SUBROUTINE FUN1 (X,FOFX> 0 1 

INTEGRAND OF STEWARDSON TRANSFORMATION INTEGRAL FOR SKIN FRICTION D 2 

COMMON AP, ZlME,E2U2,FOEtFlEK,THG12GtTWM,EMBF2t C2ME2K, X2M, WING D 3 

LOGICAL WING 0 4 

X2M=X**TWM D 5 

X?=X*X 0 6 

IF (WING) 2,1 0 7 

XF-X2M 0 8 

X2*l •/ ( X2* X2M ) 0 9 

GO TO 3 0 10 

XF-l • 0 11 

X2-X2M/X? 0 12 

F0FX-P2 ( X } **THG12G*XF*SQRT ( 1 •♦EMBF2*X2 ) 0 13 

RETURN D 14 

END D 15 


FUNCTION P (X) El 

BODY OP WING SURFACE PRESSURE 

COMMON AP,ZlME,E2U2,F0E,FlEK # THG12G,TWM,EMBF2,C2ME2K,X2M,WING*F01f E 2 

1F11K,0LR,ETB, ONETB, DUMB (131, THEM, AEP E 3 

LOGICAL WING E 4 

X?M=X**TWM E 5 

ENTRY P 7 E 6 

IF (WING) 2,1 E 7 

1 F TCE=FOF E 8 

FT1 EK=F1EK E 9 

GO TO 3 E 10 
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2 


3 


C 

A 


ETP=( ( ( DLR +ETB*XM)/(XM+( TWEM*SQRT (X)+AEP*X)*X/XM)-ET8 )*ONETB ) **2 
FTCE=F0E+ETR*(FO1-FOE) 

FT1FK*F1EK+ETR*(F11K-F1EK> 

1 me = X**Z 1 M E 

P»AP*< (U + E2U2*X21ME/S0RT(X1 >*FTOE/( X21ME+D2ME2K M-FT1EK1 


IF (PAT.l) RETURN 

PRINT 4 , P,X,E 2 U 2 , 02 ME 2 K,ETR,FT 0 E,FT 1 EK 

P *. 9999999925 

RETURN 


FORMAT (AH P -»F 13 , r fl 6 H SET * 1 * AT X =fF 10 * 5 , 10 X, 5 E 13 . 5 ) 
ENT 


E 11 
E 12 
E 13 
F IA 
E 15 
E 16 
E IT 
E 18 
E 19 
E 20 
E 21 
E 22 
F 23 - 


FUNTTION SUM(X t Y f N) 

C 

C TRAPEZ 0 TD 4 L INTEGRATION FOR UNEQUAL INTERVALS 

r 

01 KENS I ON X( 19 ),Y( 19 ) 

M=N -1 

PSUM = YMX( 2 ) -X )+Y(N)*(X(N)-X(M>) 

00 1 1 = 2 , M 

1 P SUM=P SUM+Y ( I )*<XU+L)-XC 1 - 1 ) ) 

SUK=.**P$UM 

RETURN 

END 


F 1 


F 2 
F 3 
F A 
F 5 
F 6 
F 7 
F 8 
F 9 


C 

f 

r 


1 

2 
3 


6 


7 


3 

9 


SURROUTINE SUMMARY* A, B,C ,D,E, F, H,NI 
COLLECTION OF SUMMARY RESULTS ON A FILE (TAPE 7 ) FOR SEPARATE OUTPUT 


01 PENSION 0(2), E ( 2 ) » H(8), L(2*)t 0(1A,?5) 
OAT A L,T,SKIP/26*0,3H(/)/ 

IF ( N-3 ) 1 , 2 t? 

IF (UD.EQ.il GO TO 6 
T=G 

Lim> = l 
1 = 1 + 1 
on, n=A 
0(2, n=R 
OH, I ) = c 
n (A, f ) =o 
0(5, T J=0(2I 
0(6, n=F 
0(7, n = F(2 J 
0 ( 8, 1 1 = F 
0 ( 9 ,n=H(?) 

nm,II=H( 5 ) 

n( 12 ,T)^H(A) 

0(13,1 )=H(7) 

0(lA,I )=H( B) 
tf (T.LT.25) RETURN 
FNTPY ppntsum 

WP T TF ( 7 ,9 ) 

on 7 j=! , i 

IF (L(J).PQ.l) WRITE (7, SKIP) 

WRITE (7, 6 I ( 0( K t J ) ,K®1 , IA) 

L( J) = 0 


T =0 

TF (N* 1 .TV 3 > GC TO 2 


return 

cnpM&T ( X 2 F 9 . A, F 8 . 5 , F 9 . 5 , F 8 . 2 ,F 9 . 4 ,F 8 . 2 ,F 8 # 5 , 3 X, 6 A 10 ) 

FORMAT ( 1 H 1 / 16 H HYPAEPO SUMMARY // 3 X , 5 H ( L/D ) ,*X, 5 H ALPHA, AX, 2 HCL, 3 X , 

1 2 ( 5X , 2HCD, 5 X , 5HALPHA 1 , AX, 2HCL /7 X, 3HM AX , 2X , 2 ( 3X t 5HL /OMX) , 

2 '?(6X,1H0,3X),2(2X,3HMIN,^X),3HA=0) 


ENC 


G 1 


G 2 
G 3 
G A 
G 5 
G 6 
G 7 
G 8 
G 9 
G 10 
G 11 
G 12 
G 13 
G IA 
G 15 
G 16 
G 17 
G 18 
G 19 
G 20 
a 21 
G 22 
G 23 
G 2 A 
G 25 
G 26 
G 27 
G 28 
G 29 
G 30 
G 31 
G 32 
G 33 
G 3 A 
G 35 
G 36 - 
G 37 
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Input 

A single case consists of the determination of the aerodynamic coefficients over a 
given set of angles of attack. The first card for each case provides a heading for the 
printout; it consists of 80 columns of any desired FORTRAN characters. The remaining 
cards for each case are interpreted by a system loading subroutine (NAMELIST) which 
is very flexible. The data block begins with an arbitrary name ($DATA in the present 
case) and ends with the dollar sign ($); the variables between may be in any order and 
need appear only if values are to be different from those preassigned or used in the pre- 
vious case of the same computer run. Column one of all these cards is blank. A descrip- 
tion of the input FORTRAN variables with their correct type and preassigned values (if 
any) in parentheses is as follows: 

Description 

free -stream static temperature, T^ °R (real) 

free -stream Mach number, (real) 

number of angles of attack, maximum of 11 (integer) 

angle-of-attack array, decreasing order, deg (real) 

power-law exponent, m (real) 

body fineness parameter, f (real) 

Reynolds number based on body length, R*,^ (real) 

ratio of reference area for coefficients to body 
planform area; if zero, program uses wing 
planform area (real;0.) 

ratio of x location of moment reference center to 
body length (real) 

ratio of y location of moment reference center to 
maximum body radius (real) 


FORTRAN variable 

TINF 

AMCH 

NALF 

ANGL 

EM 

F 

REL 

SSB 

XCG 

YCG 
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FORTRAN variable 


Description 


PBPI 


CAFCTR 


XOUT 


PBBPOL 


CFDCF01 


BDYONLY 


NCASE 


ratio of base pressure to free-stream static pressure 
(real;l.) 

multiplication factor times calculated laminar skin 
friction (real;l.) 

extra output at each angle of attack if 
XOUT = .TRUE. (logical, .FALSE.) 

array of NALF values of wing upper surface pressure 
parameter ^ ^ — — in ref. 2oj corresponding to 
angles of attack ANGL (real) 

array of NALF values of wing skin-friction parameter 

| — 1 in ref. 20 J corresponding to angles 

\ F,0,cr / 

of attack ANGL (real) 

set equal to .TRUE, for axisymmetric body only, 
.FALSE, for half body with wing (logical; .FALSE.) 

indicator for each additional case of a run to avoid 
unnecessary recomputations (integer; 1 initially, 

4 each case thereafter). After first case of a run 
set NCASE = 2 if AMCH or TINF is changed; set 
NCASE = 3 if EM is changed but AMCH and TINF 
are not; for no change to AMCH, TINF or EM use 
preassigned value 4. 


Output 

There are four possible output blocks for each case, only two of which always 
appear. First comes the input list with four added variables. These are GAM, the 
ratio of specific heats y; ALAMCR, a parameter (X cr ) from reference 20; DEL, the 
slender body parameter 6; and EPS, the shock strength parameter e r Next is a 
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list of the angles of attack with the pressure coefficient., Mach number and sine squared 
of the shock angle for oblique shock, and with the pressure coefficient and Mach number 
for Prandtl-Meyer expansion through the angle (appears only for new Mach numbers). 
Third is a list of variables used in the angle-of-attack and skin-friction calculations, 
which appears only if called by setting XOUT = .TRUE, in the input list. Fourth is 
the standard output of stability-axis and body-axis coefficients with the interpolated 
(L/D) max . The normal-force coefficient is also broken down into contributions from 
the body CNB, the body boundary- layer area under the wing CND, the rest of the under- 
wing area CNW, and the wing upper surface DP/Q. The axial force is broken into the 
contribution from the pressure CAP and from the skin friction CAF. In addition to these 
results, after all cases have been run, a summary of results is printed out for cases with 
angles of attack. 

Example 

Input cards for a run of two sample cases are presented below. The first case is 
for the complete configuration with m = 0.75, f = 7, at M ffl = 12, t * 256.26 x 10 6 
and 11 angles of attack. The second is for the axisymmetric body having the same 
parameters. 

POWFR-LAW TRANSPORR \ TER MACH 1 ? EM= # 7S00 F=7.00 Pr Y-4 • /m Pp Aqp -p t hjr I M T T 

SOATA AMCH=1?. •cM=.7S,F=r? # ,P(a = 256*?*F6.TINFz:40R* • Yr<i= # 1 5* 

CFDCFO 1 * G 7 * 3R * J ^ 5^ • 

NALF=1 1 * ANGL=3. 

P0WER“LAW BODY OF RpVOL MACH 1? FM=.7P00 F=7.00 ppY-4«4?F6/M PqaqF-ptmfTMJT 

$04 T A NCAS£ = ?«ROYONLY = *TPt/r . , NAt F^) , ANcL (l)T=0**PBRPOl , C F OCF0 1 =0 # * 

The output for these input cards is shown below. The total computation time on a 
CDC 6600 series computer at the Langley Research Center was less than 15 seconds 
(excluding compilation). 

Nr ASP & i, 

GAM = C.^F + Cif 

TTNF 

Nil F 
ANGl 

EM 

F 

PFL 
SSB 
x<*r, 
vcg 

PBPT 
CAFACTR 
ANG0 
XCUT 


- c.^nsF+r^, 

= O.I^F+C?, 

= Ur 

= 0.3F+01, 0.U*<n, O.lFf'M, -t. K E + OCr 0.0 f -C.5F + C0, -0.1E+01* 

-0.4. c *CJU -O.FFfO^r 

= C. 7 re + cr , f 
* 

= 0.? c '26F+n<3 t 

- r.n, 

= C.6 F*tC» 

- o. 15 f+co, 

- O.lF+0!, 

- 0.1F+0' 1 , 

■ -' 5 .U 17 f 877 lF 4 T 4 ! 5 FFC 0 t 
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PBBPOl = 0.28frF+Gi f o^lgP+nit 0.3P2e*0l f 0*?71E+01t 0.39F+01, 
O.^c + ri, 0.^?7F>0\, 0.'52F+01, 0. 4?? E+oit 0.524E+01t 
o. ee ^o», 

CFDCF01 = -0« TF-Qj p Q.6E-0W 0.23E+00, 0.32E+00t 0.42E+00, 0.52E+O0t 

Q*62 p +C0 f O.04 p +OOt 0*1 09F+O l p 0.135E+01. 0.159F+01. 

4LAMCR = C,9j96e?2l74* r 57F-C?t 

BnVPNl Y = ft 

OFI = 0. l6??52245'=4277F*O0t 

EPS = C. ?6C9ft8f 9?0f29F*PQp 

XSE * ' , -4*F-c* t 

SEND 


P0WFP-L4W TRANSPORTER MACH 12 EM=*7500 Fe7.00 FE Y*4. 42E6 ✓M PBASE-P JNFINtT 

3*00 DEG. PUPTDI = -.0060* 

2*00 OEG. PUPTOT = -.00454 

1.09 OEG. PUP 1 01 = -,002 e 7 

.*0 OFfi. PUPTOT = -.00137 

-* K 0 DFG* PUPTOT « .001*1 

-1*00 OFT,* PUPTOT s .OO?^ 

-7*00 DEC'* PUP T OT - .00746 

-3*00 O p G * PUPTOt * .0126? 

-4*00 OEG. PUPTOT = *01887 

- E .0O DFG, PUPTOT - .0263? 


(L/njWAX = 5.B*68 tT ALPHA = -.022* DFGP FF S ♦ WITH CL ANO CD = *034100 *005832 

TOO * * 0028 2901 « AT ALPHA 0 = -3.3872 

CO MIN * *00282*19, fT ALPHA a -7. WO »N0 CL =-.001197 


ALPHA 

ri 

CO 

CM 

LAP 

CN 

CN9 

CNO 

CNW 

DP/q 

CA 

CAP 

CAF 

?* 0 

.06419 

.01300 

.00211 

4.91 

.06*78 

.0*397 

.00011 

*01492 

.00578 

.00970 

.00926 

.00043 

2*0 

.0*418 

.01021 

.00024 

*.30 

.05*50 

.03764 

.00011 

.01250 

*00425 

.00832 

.00790 

.00042 

1.0 

."4424 

.00701 

-.00192 

5,6 6 

.04437 

.03174 

.00011 

*01027 

.00225 

. 00704 

.00663 

.00040 

.1 

* O'* 9 ? 8 

,00678 

-.00314 

5.79 

.03934 

.02897 

.00011 

.00923 

.00103 

.006*4 

.00604 

.00040 

C,n 

,0?4?? 

.00507 

-.00447 

5. 8* 

.03432 

,02632 

.00011 

.00825 

-.00036 

.00507 

. 00548 

.00039 


.C2°7* 

.00507 

-.00593 

5.78 

.02931 

.02381 

,00011 

.00731 

-.00192 

•00533 

.00493 

.00039 

-l.o 

.0?4?6 

.00440 

-.00754 

5. c 4 

.02428 

.02142 

.00011 

.00644 

-.00369 

.00482 

. 00444 

.00038 

-2.C 

.01428 

.00340 

-.0112* 

4.20 

.01415 

.01708 

.00011 

.0048* 

-.00788 

.00389 

.00352 

.00037 

-3.0 

.00400 

.00209 

-.0**7? 

1.36 

.00384 

.01332 

.00011 

.00347 

-.01307 

.00310 

.00273 

.00036 

-4.0 

-.00**4 

.30709 

-.02105 

-2.26 

-.00672 

.01018 

.00011 

.00234 

-.01935 

*00243 

.00207 

.00036 

*. n 4 f\ 

-.01737 

.00343 

-.02731 

-5.06 

-.01760 

.00767 

.00011 

.00144 

-.02682 

.00190 

.00155 

.00035 


SOATA 

NCASP = 2 i 

GAM = 0.14E+Clt 

TINF = 0.408E+C3t 

AMCH * 0. I 2f +02 f 

NA1 F = 1 , 

ANGL * O.Ot C.2E+P1, 0.1E+01, 0.5E+OCt 0.0, -0.5E+D0, -O-tF+Ol, 

-0.2E+01, -O.^E+Ol, -0.4E+0 1, -O.iE+Ol, 

EM = 0.7*E+CO t 

F - 0.7F+OU 

PEL - 0.2F676E+09, 

sse » O.Ot 

XCG ' = 0.6F+90, 

YCG = 0. T5 p +C0t 

PBPI = O.lE+OIt 

CAFACTF = 0.XF401, 

ANGO * -0. U17E877154745E«-O0p 

XCUT = F f 

PBBPOL * 0.0» 0.319E+01 , 0.352E+01, 0. 3718+01, 0.39F + 0W 0.409E+01, 
0.427E*01» 0.462F*01t 0*49*F+01t 0-524F+01, D.55F+01t 

CFDOFOl = O.Ot C.6F-01, 0.23F+00, 0.32E + 00, 0.42F+C0t 0.52E+O0, 

O.A2F+OOt O.84F+00, 0.109F+0l t 0.135F*0l t O.l^SE+Ol, 
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APPENDIX - Continued 


ALA MCA * 0.9196S2m*6557E-C2„ 

6DYONLY * T, 

OEL = 0.1632522A594277E+00, 

EPS - '1.26O*66e ,? 920f?5E + Ee* 

XSE * 0.1?«3960a*l*ll69F-0? t 

SEND 



powfr-lah pcdy 

OF REVOL 

MACH 12 

EMa.79C0 

F*7.O0 

PPY=***2E6/M 

P0ASE 

-PINFINJT 




ALPHf 

CL CO 


L/D 

CN 

CNB 

CND 

CNW 

0P/0 

CA 

CAP 

CAF 

0.0 

0.0W0 .01128 0. 

00000 

0.00 

0.00000 

.07632 

.OOOIi 

.00025 

0.00000 

.01128 

.01096 

.00032 


HVPAfPO 51 WM 4 HY 


<L/P) ALPHA rt CO 

MAX l./DMX t/DMX 0 


ALPHA CO ALPHA CL 

0 MTN pin a*o 


5.8*68 • -.0?2* : .03A\0 .00283 -3.38 


.0028 -3.50 .03*32 TEP MACH 12 «=M«.75aO fsT.OO RFY**.*2E6/M P8AS6=P I NF INI 7 


COMPUTER PROGRAM FOR CALCULATING Kg AND X cr 


The main program requires as inputs values of two parameters describing wing 
upper su 2 'face conditions. These values are an average wing upper surface pressure 
parameter (PBBPOL) and an average wing upper surface skin-friction parameter 
(CFDCF01)._ They are plotted in reference 20 (figs, 4 and 11 of the reference, respec- 
tively) as — - — ^ and F,A — 1 for delta wings as functions of (a viscous 

X cr C F,0,cr Cr 

interaction parameter) and K Q ( = -M w a). This program calculates K 0 for each angle 
of attack and the value of A cr for the delta wing corresponding to the power -law wing. 
As mentioned in the main body of the paper, the correspondence is based on the viscous 
effects, which are assumed to be approximately equal for wings with equal spans and 
equal values of x dr dx. For the power -law wings, this integral involves gamma 
functions which are approximated analytically in the program. 


l 

? 


PROGRAM UPPPESSC INPUT. 201, OUTPUT* 201, TA PEI = INPUT I 
DIMENSION ANGL(ll), C«0U1), P0(11J, HE ADI 8 ) , P80POU11) 

NA«ELTST /DATA/ RE1 , F , AmCH , T t NF , EM , FT A6,F0,A1,NCASE, NALF , ANGL, X r G* 
ISSB, Y<“G 

DATA GAM ,TWTT,PTf, ASAV/1. A,. 41687,2.68868,0./ 

GM1.GAM-V . 

GMlf.l./GMl 

GM83=GAM*8./3. 

GM12=.**GM1 
GPI 4=1 GAM*1« ) /4. 

ZGGI=GAM/GM1 2 
READ 11, HEA0 
IF (ENDFILE 1) 9,2 
READ OATA 

IF ( ANGl.EQ.ASAVI GO TO 7 

ASAV.ANGL 

DO 8 1 = 1 , NAL F 

ANG=. 0174533* ANGtl I I 

CAOI I )=-AMCH*ANG 


A 1 
A ? 
A 3 
A 4 
A 5 
A 4 
A 7 
A 8 
A 9 
A 10 
A 11 
A 12 
A 13 
A 14 
A 15 
A 16 
A 17 
A 18 
A 19 
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APPENDIX - Continued 


3 


L 


6 

7 


8 


9 

r 

10 


n 

12 


IF <CA0U)) 4,5,3 
GPlK04=GP14*CA0<I1 
PTGK4*SQPTf i .*GPlK04**2) 

POH l*l.+G4M*C40< T)*(GP1K04+RTGK4) 

P3BP01_( T)=GM03*< RTGK4+ GP IK 04* { 2 • +GP1 K 94/P TGK4 ) )/5QPT( PO (1)1 
GO TO 6 

GM12K*l.+GM12*CA0n , 

POtT 2K** ZGG1 

IF (P0( n.LT.O.) P0<T)*O. 

PB0P 01 ( T)=GM83*GM12K**GM1T 

GO TO 6 

pon i«i # 

P86P01 < T > =GM83 

CONTINUE 

ZM=?.*pm 

PRINT 12, HEAD 

IF (FM.FQ.B64V) GO TO 8 

BSAV-F* 

FM4019=„4019*EM 
FMXl^l . r 303-FM&01 9 
FMX2=.5/( .9274+EM4019 ) 

AT*.977 B ?-* i 7627*EMM1/<EMMH>l,295*FM4l ) 

FMG-GMi*AT*ETAB/(EM*SQRT( 2.*GM*EM^i+F0) ) 

AFP*A1*(FT4B*F/AMCH)**2 
Aepx=AF0*EMXl/(5.-ZMl 
TTBT7-1 .+GM1 2*AMrH*AMCH 
TWTI=TWTT*TTBTT 

TPTI*.273M.19f>+.532*TWTT)*TTBTI 

SUTMT=198,6/T1NF 

Cl = ( 1. ♦ SUTHT I/(TWTI^S1JTHI )*5QRT(TWTT ) 

C3 = n .+SUTHI )/<TPTUSUTHT>*SqpT(TPTT ) 

ZEM32=FMG*4MCH*F*F*$QPT< Cl/* 5 El ) 

EM?2X=.5*ZEM*2/EMX! 

G2MPC=.*90F3*GMI2*<AMCH**3 ) *SQPT < C3/R El) 
rpE!.,=PTF*<FMX2 + 4EPX+FM32X) /<1. *AEP + ZEM32) 

ALAMCP«g?MRC /SQPT fCREL J 

PRINT IOt EM,F, AMCW»C°EL, ALAMCR, ( ANGL ( T >,CAOU ),POI T ) ,PBRPOL(I), 1= 
lit NAIF) 

GO TO 1 
STOP 

FORMAT t///15X f 2HFM T 9X f 1.HF f 8X,&HMArH T 6X t FHCPE/L f 7X t llH{L AMPOA)CPF/ 
l/9X,Fin*4,2F!O.2,F12.6,ElT*6///6X > r H'\NGlF,6X,4H(K)0,7X,AHmO,?X f 6 
2HI AM=0. ,9X, 13H l P-PO) /LAMBDA, 6X,1'H( CF C/CFO 1-1/1 /FI 1.2 , FI 1.4, JFll.* 5 
3>> 

FORMAT ( 8 A 1 0 I 

FORMAT C1H1» 5X8410) 

END 


A 2 0 
A 21 
A 22 
A 23 
A 24 
A 25 
A 28 
A 27 
A 28 
A 29 
A 30 
A 31 
A 32 
A 33 
A 34 
A 35 
A 36 
A 37 
A 38 
A 39 
A 40 
A 41 
A 4? 
A 43 
A 4* 
A 45 
A 46 
A 47 
A 48 
A 49 
A 50 
A 51 
A 52 
A 5? 

A *4 
A ^5 
A 56 
A 57 
A 58 
A 59 
* 60 
A 61 
A 62 
A 63 
A 64 
A 6 P 
A 66 
A 6^ 

A 68 
A 69- 


Input 

A single case consists of calculations for a single configuration over a set of angles 
of attack* The first card is a heading consisting of any 80 FORTRAN characters* The 
remaining cards use the same loading subroutine as does the main program; the data 
block begins with $DATA and ends with $. The necessary input variables are TINF, 
AMCH, NANG, ANGL, EM, F, REL, ETAB, FO, and Al. Of these the first seven are the 
same as for the main program and the last three are from the similarity solution results. 
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APPENDIX - Concluded 


See figure 3 and table I where ETAB = 77^, FO = Fq(? 7 ^) } and A1 = aj. The program is 
set up so that the input cards to the main program may be used with these three variables 
added. Here is the input for the first sample case given with the main program: 


POWFR-LAW TRANSPORBlTtP MACH 1? EM=.7^00 F=7.00 PrY =4.A?F4 /m 

SOATA AMch= 1?.,EM=.7S,F = 7., T INF=40R . « XCO=.ft , YCn=. 1 S » 

ETAB=.87S07,F0 = .(S9A0ft.Al = ,ft071?, 

NALF = 1 1 ♦ ANfiL*3. ♦ 1 • ♦ »5«0« * ,-o, »-3. $ 


« ^ * flPp^F^c; 


Output 

The output for the same case is shown in this section. In it (LAMBDA)CRE = x cr 
and (K)0 = Kg are needed for use with reference 20. Also printed are the length ratio 
of the delta to the power-law wing, CRE/L; the ratio of inviscid surface pressure to 
free-stream static pressure, PO; and the average pressure parameter for X = 0, 

LAM = 0. This latter value is a useful aid sometimes in interpolating values from the 
figures of reference 20. 


POWFR-LAW TPANSPOORITFR MSCH 12 FM=„7*00 F*7.00 PEY=4.42Ff/M ***** (JPPPFSS 



EM 

F 

MACH C9F/1 

l LAMBDA JCRf 


♦ ^600 

7.00 

12.00 1.056144 

9.196840?-')? 

ANGI.F 

<K JO 

< p 1 0 

t AM = 0 . 

<P-P0) /LAMBDA 

3.00 

-.628^ 

.39062 

2.66866 


2.00 

-.41 89 

.54202 

2.99986 


1 *00 

-.3094 

.74116 

3.35458 


.*0 

-.1047 

. 862?9 

7.54092 


0.00 

-0.09OQ 

1.00000 

3.73333 


-.50 

.10*7 

1.15611 

3.92899 


-1 .00 

• 209* 

1.33737 

4.12332 


-2.00 

.4189 

1.75205 

4.49871 


-3.00 

. 62 83 

2.27170 

4.84A14 


—*•00 

.83^8 

2.90224 

5.15051 


-5.00 

1.^472 

3.65262 

5.41*72 



I C F0/CF0) -1 
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TABLE I.- SOLUTIONS TO THE HYPERSONIC SIMILARITY EQUATIONS 
FOR POWER-LAW BODIES OF REVOLUTION 


m 

% 

F o(%) 

F i M 

J 0 

J i 

a i 

r - 7/5 J 

1.00000 

0.91492 

0.87342 

0.9179 

0.07323 

0.08410 

0.47546 

.95000 

.91034 

.84711 

1.0591 

.07589 

.09551 

.52709 

.90000 

.90465 

.81630 

1.2306 

.07909 

.11044 

.58604 

.85000 

.89743 

.78174 

1.4386 

.08303 

.13067 

.65291 

.80000 

.88798 

.74265 

1.6887 

.08799 

. .15918 

.72741 

.75000 

.87507 

.69806 

1.9811 

.09444 

.20098 

.80732 

.70000 

.85648 

.64662 

2.2986 

.10318 

.26459 

.88631 

.66667 

.83880 

.60763 

2.4964 

.11098 

.32565 

.93216 

.63333 

.81391 

.56403 

2.6392 

.12129 

.40794 

.96566 

.60000 

.77647 

.51478 

: 2.6593 

.13564 

.51763 

.98034 

.55000 

.66414 

.42678 

2.2510 

.17318 

.74598 

.96791 

.53000 

.56901 

.38500 

1.8876 

.20119 

.86687 

.96377 

.51000 

.37221 

.33757 

1.4110 

.25443 

1.04110 

.97539 

.50500 

.27299 

.32450 

1.2766 

.28069 

1.11845 

.98249 

.50000 

.00000 

.31077 

1.1366 

.35808 

1.36841 

.99182 

y - 5/3 

1.00000 

0.87041 

0.81065 

0.7836 

0.10244 

0.10987 

0.46531 

.95000 

.86429 

.78363 

.9017 

.10532 

.12597 

.51356 

.90000 

.85679 

.75282 

1.0433 

.10872 

.14660 

.56788 

.85000 

.84740 

.71823 

1.2122 

.11283 

.17364 

.62833 

.80000 

.83532 

.67912 

1.4108 

.11787 

.20994 

.69402 

.75000 

.81919 

.63448 

1.6356 

.12422 

.25974 

.76228 

.70000 

.79658 

.58296 

1.8685 

.13248 

.32919 

.82727 

.66667 

.77569 

.54389 

2.0053 

.13956 

.39031 

.86398 

.63333 

.74719 

.50016 

2.0956 

.14850 

.46626 

.89116 

.60000 

.70595 

.45067 

2.0942 

.16025 

.55904 

.90627 

.55000 

.59076 

.36177 

1,7872 

.18792 

.73266 

.91473 

.53000 

.49985 

.31912 

1.5217 

.20647 

.81845 

.92427 

.51000 

.32217 

.26988 

1.1590 

.23916 

.93877 

.94825 

.50500 

.23542 

.25600 

1.0506 

.25495 

.99164 

.95764 

.50000 

.00000 

.24113 

.9315 

.30378 

1.16431 

.96872 
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(a) General power -law, wing -body configuration. 



(b) Body of revolution in nondimensional coordinates. 

Figure 1.- Configuration studied, showing relation between physical 
and nondimensional coordinates. 
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Slender body parameter, d 


Figure 2.- Graph of the relation e. = l/(M oo 6)^ for several Mach numbers and 
showing the values of 6 and e-, for various sets of experimental data. 





in 

HI SB 

I BHI 

■ \mi 




Figure 4.- Effective body for estimating aerodynamic characteristics at 
angle of attack (shown for f = 6, m = 2/3). 
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Theoretical drag forces 

* Total drag (present method) 

Hypersonic similarity 

Skin friction 

Newtonian method 

t n v i sc id cone flow 



Power-law exponent, m Power-law exponent, m 

Figure 5.- Comparisons of theoretical drag coefficient at zero angle of attack with 
experiment for power-law bodies of revolution. The theory is shown for the 
larger Reynolds number in each case. 
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Angle of attack, a, deg 


Figure 6.- Comparison of theoretical normal and axial forces with experiment for 
wing— conical -body configurations {m = 1) at Mach 10.03. 





Figure 7.- Comparison of theoretical drag polar, lift-drag ratio, and pitching- moment 
curves with experiment for wing— conical -body configurations (m = 1) at Mach 10.03 





(a) Drag polars and pitching moment at Mach 6.86 in air. 


Experiment (ref. 7) 



(b) Maximum lift-drag ratio at Mach 19.9 in helium. 

Figure 8.- Comparisons of calculated performance of half-cone bodies having 
shock- shape -matching wings with experimental data for the same bodies 
having delta -planform wings. 
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